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ABSTRACT 


The expression for the exospheric temperature in Jacchia's 
static diffusion models of the upper atmosphere has a discon- 
tinuous gradient at the poles. Therefore it cannot describe the 
true state of the upper atmosphere in the polar regions. Further- 
more, it cannot be used to calculate quantities that depend on the 
derivative of the exospheric temperature, or the density, like 
pressure gradients, horizontal forces or horizontal heat fluxes. 

A modified expression for the exospheric temperature is sug- 
gested. This modification yields variables of state of the upper 
atmosphere that deviate little f:om Jacchia's values, but it has 
continuous gradient at the poles and is therefore more suitable 
for treating dynamical problems like the global wind pattern. 


ON EMPIRICAL MODELS OF THE UPPER ATMOSPHERE 
IN THE POLAR REGIONS 


INTRODUCTION 

In the past decade and a half a very great number of observations of the 
physical state of the upper atmosphere were obtained both by satellite drag 
analysis and by direct measurements with rockets and satellites. The analysis 
of these data led to the development of Jacchia's semi-empirical static diffusion 
models of the upper atmosphere (Jacchia, 1965 and 1971). These models de- 
scribe the observed variations of the upper atmosphere generally with great 
accuracy. In Jacchia's description of the upper atmosphere all state variables 
depend upon the exospheric temperature, i.e., once the exospheric temperature 
is determined, the state of the upper atmosphere at any height can be ascertained. 

According to Jacchia the exospheric temperature depends upon the solar 
activity F, the average solar activity F during several solar rotation periods, 
the geomagnetic activity index K p , the day in the year (semi-annual effect), the 
local solar time (diurnal variation) and the latitude. The dependence on LST 
(local solar time) can also be expressed as an azimuthal dependence, where the 
azimuthal angle <p is measured from the azimuth that corresponds to a fixed 
LST. 


The variation of exospheric temperature with local time and with latitude 
can be given by a single mathematical expression, i.e., the distribution of exo- 
spheric temperature on the globe, where the longitude is measured with respect 
to local noon. 

In Jacchia's model the exospheric temperature is given (in his notation) by: 


T c (1 + R sin m 0) 


( cos m v - sin m <9 

1 + R 

\ 1 + R s in m Q 



where T c depends only upon solar activity and 


T - H + /3 + p sin (H + y) (- n < t < 77) 

_ \<t> ~ Si 
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0 - latitude 

S - declination of the Sun 
II - LST counted from local noon 
(i = -37° 
p = G° 
m = 2.2 
R = 0.3 
y = 43° 

For the problem treated in tins paper we shall use the notation: 


T a " Fj(£) + F 2 (tf) • F 3 (<p) 


(2) 


where i) and cp are the colatitude and the longitude respectively, and 


F,(fl) = T C (1 + R sin m i9) (3) 

f 2 (i>) r T c R(cos m 7j - sin m (9) (4) 

F 3 (<P) = cos3 (^J (5) 


The subject of this paper is to show that expression (1) for the exospheric 
temperature cannot describe the physical conditions near the poles in a true 
fashion. We shall suggest a modification of (1) that deviates very little from it 
at all latitudes, but wiL give more consistent results in the polar regions. 
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JACCHIA'S EXOSPHERIC TEMPERATURE 
FOR THE POLAR REGIONS 

Jacchia's model is a static model; this implies that it is a description of the 
steady state of the atmosphere. The global variations of temperature and density 
that result from it should therefore be continuous functions of the independent 
variables like height, latitude and local time. In nature there appear no discon- 
tinuities of the variables of state or their derivatives in the steady state. As in 
Jacchia's model all variables of state depend upon the exospheric temperature, 
a discontinuity of it, or of its derivative, will cause also discontinuities of the 
other variables or their derivatives. We shall show in the appendix that expres- 
sion (1) has a discontinuous derivative at the poles and therefore cannot be a very 
good approximation to the state of the upper atmosphere in the polar regions, 
although a direct comparison of Jacchia's model with observations in the polar 
regions is difficult, as very few exact observations exist for these regions. 

By a derivative of a variable defined on a sphere we understand the deriva- 
tive of this variable along any arc of this sphere. The value of this derivative 
depends on the particular arc cnosen. By a continuous derivative of the variable 
at a point of the sphere we mean that the derivative of the variable along an 
arbitrary arc is continuous. 

The discontinuity of the derivative of the expression for the exospheric 
temperature at the poles can easily be demonstrated by formal differentiation. 

In this paper we will not do this, but shall .analyse (in the appendix) expression (1) 
in terms of spherical harmonic functions, as our suggested modification becomes 
more transparent by this method. 

All single valued functions that are defined for all points of a sphere given 
by colatitude d and azimuth cp(0£d<77, 0 £ cp < 2??) that are continuous and 
have continuous derivatives can be developed into a series of spherical harmonic 
functions. This series will converge uniformly for all points of the sphere to the 
value of the function and may be written 

F(d.cp) = £ £ (a mn p n n ’(cos d)cosmcp + b m[ P n m (cos fl) sin mcp) (5) 

n=0 m=0 

Jacchia's function cannot be developed into such a series for all points of 
the sphere. It has a discontinuous derivative at the pole, or what is equivalent, 
a cusp at the poles. This cusp can easily be removed by modifying the term 
cos 3 (t/2) that appears in the original expression. By making the difference 
between the modified term and the original term small, the resulting exospheric 
temperature will deviate only little from Jacchia's original values. 
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In order to reach this we suggest to substitute for cos 3 (t/2) the expression 


o • 2 q /l + cos (cp - 37 )\ , 3 / x \ 

S(i3. <p) = sin 2 d • I 2 j + cos 2 d cos 3 


(6) 


At the equator both expressions are exactly identical. At the poles the function 
F 2 (d) is equal to zero, so (6) does not influence the value of the exospheric 
temperature at all. On the other hand the derivatives of the exospheric tem- 
perature become continuous by this substitution. 

Our suggested expression for the exospheric temperature will be in Jacchia's 
notation 

( cos"’?? - sin m (9) 

1 + R - 

1 + R sin m £ 

( , ( 1 + cos ( 0 ) - 37°) \ , / t \\ 

I COS 2 4> ^ 2 ) + Sin 2 ^ cos 3 ^ J J 

or in our notation 


= T .(l + R sin m $) 


T t Fj(i^) + F 2 (0)S(d.q>) 


The physical meaning of our modification compared to Jacchia's original 
expression is to leave the diurnal variation unchanged at the equator. We retain 
the first order Fourier coefficients at all latitudes, but reduce the zero and the 
higher order Fourier coefficients gradually to zero as the latitude approaches 
the pole. 

Furthermore, the substitution will not cause any shift in the position of the 
global maxima and minima of the exospheric temperature, i.e., the modified 
function will retain the property that the maximum of the temperature coincides 
in latitude with the subsolar point and shows a phase shift of 37° in longitude. This 
is easily seen as at the extrema cos 3 ( T . 2) becomes either zero or one and 
cos (q? - 37 c ) either 1 or -1. Substituting these values in (6), the expression be- 
comes independent of latitude and therefore does not influence the derivative of 
T, with respect to latitude at the longitude of global extrema. The position of 
the extrema is therefore unchanged. 
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Tables 1 and 2 show the ratio of the modified temperatures to Jacchia's 
original values. Nowhere is the deviation greater than 1.7%. The exact shape 
of the diurnal variation is only changed very slightly, less than can be verified 
by observations. 

Figures 1 and 2 show the lines of equal exospheric temperatures in the polar 
region. It is immediately apparent that the line that passes through the pole ac- 
cording to the original model has a sharp edge there. 


ATMOSPHERIC FORCES 

The horizontal driving force in the atmosphere is given by 

F h = - (grad h p) / density, 

where grad h is the horizontal gradient. 

We have calculated the driving forces according to Jacchia's model and 
according to our modification. In the first case we have found that the driving 
force jumps at the poles by a factor of 3.5 and at some longitude even shows a 
change in sign. Our model results in smooth forces at the poles. 

In calculating the global wind field it is essential to use a well defined smooth 
force everywhere on the globe, as discontinuous forces at the poles cause the 
calculated results to be distorted at all latitude, because the boundary condition 
can influence the values in the whole domain. 


RESULTS 

(1) Jacchia's semi-empiric model of the upper atmosphere does not give a con- 
sistent description of the state of the atmosphere in the polar region. 

(2) Especially, the model is unsuitable to calculate pressure forces, heat fluxes 
and gradients of variables of state in the polar region. 

(3) A modification of Jacchia's formula for the exospheric temperature distribu- 
tion will make the model smooth in the polar regions and will deviate only 
slightly from the exospheric temperatures given in Jacchia's original ex- 
pression at any latitude. 
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(4) The horizontal pressure forces derived from the modified model deviate 
considerably from the forces derived from Jacchia's model. At high lati- 
tudes the difference in amplitude can amount to 30% of the global maximum 
force. 
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APPENDIX 


A single valued function F(d, <p) that is continuous and has continuous deriva- 
tives at all points of a sphere can be developed into a uniformly convergent series 
of spherical harmonic functions: 


F(d.<p) = £ £ I a mn P n m (cos d) cos(mcp) + b mn P n n '(cos d) sin(nnp) 

n = 0 m = n L 


(la) 


or as re-arranging terms of a uniformly convergent series is allowed: 


F(d,<J>) = £ 


ra = 0 


£ . 

m = 0 L 


£ a mn P ,r( cos ^) cos ( mc P) + b mn P n m (cos d) sin (mcp) 


A m (cosd) cos (mcp) + B m (cosd) sin (mcp) 


(2a) 


where B 0 (cos d ) is defined as identical zero. 
We shall use the notation 


p ;' <“» » > = p " (cos » > 


(3a) 


The spherical harmonic functions have at the poles the following properties: 


P° (pole) = (±l) n ; P" (pole) = 0 if m * 0 

P m ' (pole) = 0 if m ^ 1 


(4a) 


(5a) 


therefore: 


A 0 (pole) may have any value, 


(5a) 
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(6a) 


A ra (P°le) = B„,(Po1p) 
A m ' (pole) = Bj (pole) 


0, if m \ 0 , 
0, i f m t 1 , 


Aj (pole) and Bj (pole) may have any value. 


(7a) 

(8a) 


As examples of functions on the sphere that have discontinuous derivatives 
we may take sin d and sin d cos 2 cp. For the first function A Q = sin d and 
therefore A 0 ‘ = cos 0 which is not equal to zero. In the second case the same 
argument applies to A 2 . This shows that a development of these functions into 
a uniformly convergent series of spherical harmonics is not possible, they can- 
not belong therefore to the class of functions that are continuous and have con- 
tinuous derivatives for all points of the sphere. In fact, these functions have a 
cusp at the poles. 

We shall show that also Jacchi i’s function is of this type, it is continuous at 
the poles, but has a cusp there. 

In our notation 


T = F,(d) + F a (d) * F 3 (cp) (9a) 

By a numerical analysis of F 3 (cp) = cos 3 ( T /2; we obtain: 

F 3 (cp) = 0.4194 + 0.4088 cos(cp) + 0.2984 sin(cp) + U(<p) (10a) 

where U includes all the higher terms. Their coefficients are decreasing to zero 
like n -4 . The first order terms may also be combined 

0.4088 cos (cp) + 0.2984 sin (cp) = 0.5061 cos (cp -37°) (Ha-) 


Therefore 


Fj(fl) + F 2 («*) ' 0.4194 + F 2 ($) • 0.5061 cos (cp -37°) 


(12a) 


♦ F 2 (tf) • U(cp) 
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We can now identify: 


A 0 (cosd) 

F,(d) x 0.4194 

Aj (cos d ) = 

0.4088 F 2 (d) 

Bj (cos d ) = 

0.2984 F 2 (d) 


(13a) 


while the A . and B . with i greater than 1 have the same d -dependence as A r 
At the poles they behave like F 2 , which means that they are zero, but have a 
derivative that is not zero. 

Simple differentiation of A 0 with respect to d gives the result that its 
derivative is unequal to zero at the poles. The same is true for all the other 
A t and Bj. 

The cusp at the poles of the function describing the exospheric temperature 
can be easily removed if we modify the term cos 3 (r/ 2) in such a way that the 
values of the resulting A . and B . will conform at the poles to the required con- 
ditions, but the modified term will deviate as little as possible from the original 
term. 

The modification by which this is achieved is giver in the main text. 
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RATIO OF MnniFIFD EXOSPHERIC T F f'PERATURE TO JACCHIA'S EXOSPHERIC TEMPERATURE 
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LINES OF EQUAL EXOSPHERIC TEMPERATURES 
FOR POLAR REGION ACCORDING TO JACCH1A MODEL 


Figure 1. Lines of Equal Exospheric Temperatures in the Polar Region According to Jacchia's 
Model for the Global Temperature Distribution 
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